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The DELPHI and ALEPH experiments have released new results of the Fermi-Dirac correla- 
tions for pp (pp) and Bose-Einstein correlations for K^R®. Both experiments measure very 
small source radius for protons R ~ 0.1 fm. The source dimension for K®K° is in agreement 
with the previous measurements for K + K~ , K°Kg. These results, together with earlier LEP 
measurements, establish the dependence of the correlation radius on the hadron mass. This 
paper discusses some of the attempts to describe this phenomenon. 



1 Introduction 



The Fermi-Dirac (FDC) or Bose-Einstein correlations (BEC) are examined by measuring a 
two-particle correlation function denned as a ratio of the inclusive density distribution for two 
particles and the so-called reference density. The reference density is a two-particle density 
distribution which approximates the distribution without the FDC or BEC effect. The central 
problem in the FDC or BEC analyses is the definition of the reference sample. There are 
three different ways of defining the reference sample: (a) the Monte Carlo sample without the 
FDC or BEC, (b) an event-mixed sample based on the choice of two identical hadrons, each 
originating from a different data event, (c) the distribution for opposite charged particle pairs 
(the so-called "unlike" distribution). In the reported analyses the correlations are studied in 
the four-momentum difference Q, where the correlation function is fitted using a Goldhaber-like 
parametrisation^ 

C(Q) = N(l±Xe~ Q2R2 ), (1) 

where R is the source radius, A is the so-called chaocity parameter which characterises the 
strength of the correlations, N denotes the normalisation factor. In the above equation the sign 
" + " is assigned to the BEC effect, " - " to the FDC effect. 



2 Fermi-Dirac correlations for protons (antiprotons) 



2. 1 DELPHI analysis for antiprotons 2 



The data sample of ~ 2 M of Z° —> hadrons events, recorded in 1994 and 1995 in the DELPHI 
detector at LEP was used. The Monte Carlo simulation sample was obtained using the JETSET 
7.4 generator (without the FDC effect) and contained ~ 6 M hadronic events. Since it was cru- 
cial to identify the particles the requirement for the minimal particle momentum of 0.7 GeV was 
used. In this momentum region the MC agrees very well with the data and the identification 
information is available from both the RICH detector and the dE/dx measurement. After ap- 
plying the selection cuts there remained about 180 K events with at least two proton candidates 
for the analysis in data and about 560 K in the MC sample. The single proton purity was 85% 
at the efficiency of 70%, as determined from the MC simulation. This resulted in 70% purity for 
the pp sample and 50% purity for the pp sample. The cleaner pp sample is used for the source 
radius measurement. The analysis of the pp sample, after background subtraction, yielded con- 
sistent results. The correlation function was determined using three different reference samples: 
the Monte Carlo, event-mixed and unlike sign pairs. The correlation function, together with fit 
results based on parametrisation Q for pp, constructed with respect to three different reference 
samples, are shown in Fig^ Since the analysis relying on the event-mixing reference sample 
was more precise than the MC or unlike analyses, the final result was taken from the event- 
mixing analysis: R = 0.16 ± 0.04(stot) ± Om(syst) fm, A = 0.67lg;J?(stat) ± 0.18(syst). The 
other two methods were used to estimate the systematic error. The main contributions to the 
systematic error come from the proton identification, variations in different reference samples, 
parametrisation of the correlation function and binning width in Q. 
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Figure 1: Correlation function for pp (DELPHI preliminary results): (a) event-mixed reference sample, (b) MC 
reference sample, (c) 'unlike' sample. Solid lines denote fit results using parametrisation Q. 



2.2 ALEPH analysis for protons (antiprotons) 3 



The analysed data sample consisted of a total of 3.9 M Z° — ► hadrons events and 6.5 M 
Monte Carlo events obtained employing the JETSET 7.4 generator (without the Fermi-Dirac 
correlation effect). After applying the selection cuts the purity for a single proton was 70%, 
which gave a 74% purity for (anti) proton- (anti) proton pairs for Q < 5.0 GeV. The fits to the 
correlation function are shown in Figj2j4eft). The discrepancy between the proton momentum 
distributions in the data and the MC simulation in the region 0.5 < p < 1.3 GeV was corrected 
by using momentum dependent weights. The region 1.3 < p < 3.0 GeV was excluded from 
the analysis since the identification was based only on dE/dx. The analysis was conducted 
for the event-mixed and the Monte Carlo reference samples for both pp and pp pairs. The 
final result quoted for the event-mixing analysis is: R = 0.10 =L 0.02(stat) ± 0.02(syst) fm, 
A = 0.46 ± 0.04(stat) ± O.ll(syst). The main contributions to the systematic error are due to: 
different reference samples, the Monte Carlo reweighing, the Coulomb corrections. 
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Figure 2: ALEPH results (preliminary): (Left) Ci(ff) for pp and pp pairs using different reference samples: (A) 
Monte Carlo, (B) mixed event double ratio. The solid lines indicate the results of the fits using the parametrisation 
JTJ. (Right) C 2 (Q) for K°K° pairs using different reference samples: (A) Monte Carlo, (B) mixed event double 
ratio. The curves represent the results of the fits using parametrisation Q (solid lines) and the exponential 
parametrisation (dashed lines). The region of /o(1710) (open circles) has been excluded from the fits. 



3 ALEPH analysis of Bose-Einstein correlations for R^K® 

The analysis was based on the sample of ~ 200 K events with two K®. The purity of the selected 
K®K® sample was 96%, the efficiency 27%. The correlation function was obtained using two 
reference samples: the Monte Carlo and event-mixed. The fits to the correlation function 
were performed with two different parametrisations: the Goldhaber-like (pQ) and the exponential 
(FigEfright)). The long range correlations had to be taken into account for the Monte Carlo 
reference sample in the region of Q > 0.8 GeV. The presence of the /o(980) resonance, simulated 
in the MC with the width of 0.1 GeV, affects the first bin. The Q region corresponding to the 
/o(1710) resonance was excluded from the fit. The above effects account for the systematic 
error. The final result comes from the event-mixing analysis: R = 0.57 ± 0.04(stat) ± 0.06(syst) 
fm, A = 0.63 ± 0.06(stat) ± 0.06(syst). 



4 Correlation radius dependence on hadron mass 



The above results, together with the earlier LEP results^, establish a strong decrease of R with 
the increase in hadron mass. There are three approaches to understand this effect: semi-classical 
derivations based on the Heisenberg relations or the virial theorem^ and a quantum-mechanical 
model based on the Bi or ken- Gottfried condition^. 

The first attempt - links the R(m) behaviour to the Heisenberg uncertainty principle. Using 
the Heisenberg relations in momentum-space and energy-time, and identifying At with the 
time scale of strong interactions (At = 10 -24 s, independent of the hadron mass), one obtains 
R(m) = °^ 4 j , which correctly reproduces the data. 

In the second attempt^!, the virial theorem is used to estimate the relative momentum of two 
particles bound by the potential V(R). For the general QCD potential V(R) = kR — g 2 ^ 2 with 
the parameters set to the values obtained from the hadron wave functions and decay constants, 
the radius of the particle emission region changes according to the formula: R ~ 1/ ^fm. 

Both explanations suggest that hadrons are emitted not from a unique source, but from the 
sources with radii strongly dependent on the mass. This makes the energy density of the source 
emitting protons or A°'s unacceptable^. 

In the third attempt to understand R(mP, the radius of the particle source is universal. The 
observed R(m) dependence is solely the consequence of the Bjorken-Gottfried relation, which 
connects the space-time position x^ of an emitted hadron to its four-momentum q^: q ll = Xx^. 
This relation is used in a quantum-mechanical formulation of the model by postulating that the 
source function is factorised: S(P,X) = F(t)S\\S±, with the S|| and S± having Gaussian forms 
representing space-momentum correlations. The source function is related to the density matrix 
in momentum space by the Fourier transform p(q, q') = J d 4 Xe l< ^ x S(P, X) , where P = (q+q') /2, 
Q = q — q' and q, q' are four-momenta of a particle. Therefore, if the source function is known, 
the one- and two-particle distributions are determined completely. The parameters of the model: 
the ranges of position-momentum correlations, were constrained by the fit of the predicted pt 
distribution to the data 121 . The model reproduces correctly the observed R(m) dependence and 
the longitudinal source elongation. An important consequence of the model is the fact that the 
observed source radius is the apparent radius, which is smaller than the real particle emission 
radius: 

R^ s * & + # Robs < R > (2) 

where tq is the longitudinal proper time, 5 parametrises the correlation range. 
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